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METHOD FOR CALCULATING CONVECTIVE HEAT -TRANSFER 
COEFFICIENTS OVER TURBINE VANE SURFACES 


by Daniel J. Gauntner and James Sucec* 

Lewis Research Center 

SUMMARY 

A method for calculating laminar, transitional, and turbulent convective heat- 
transfer coefficients over the surfaces of turbine vanes is described. Results of calcu- 
lations with an approximate integral solution method are compared with results from a 
turbulent flat plate solution and a finite-difference solution. 

The approximate integral solution method produced heat-transfer coefficients in 
good agreement with the finite -difference solution method in the laminar and fully turbu- 
lent flow regimes. The transitional flow regime length was approximately the same al- 
though the regime began earlier for the finite -difference solution. The turbulent flat 
plate solution produced a conservative heat-transfer distribution which exceeded the 
finite -difference solution by as much as 150 percent. 

The calculations were made for given gas stream conditions of turbine inlet tem- 
perature equal to 1273 K, turbine inlet pressure of 45. 7 newtons per square centimeter, 
a fuel/air ratio of 0. 0, a vane channel inlet Mach number of 0. 31, and an exit Mach 
number of 0. 84. Transition criteria were based on momentum thickness Reynolds num- 
bers in the integral solution. 


INTRODUCTION 

In the design of blades and vanes for gas turbine engines, the convective heat- 
transfer rate from the gas stream to the surface must be calculated. The simplest ap- 
proach is to use the flat plate solution. The flat plate solution, however, applies only to 
zero pressure gradient flow. In many instances, though, this approximation leads to 
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satisfactory results for initial calculations because the Stanton number is relatively in- 
sensitive to pressure gradient (ref 1). One of the most powerful methods of calculatii^ 
the heat flux to a turbine blade is to solve the boundary layer mass, momentum, energy, 
and turbulent energy equations simultaneously by a finite-difference approach (ref. 2). 

In between these two approaches is the approximate integral solution. This method is 
more general than the flat plate solution, but not as powerful as the finite -difference 
methods. 

An approximate integral solution (ref. 3) for the local surface heat-transfer coeffi- 
cient that takes into account the effect of pressure gradient and of variable surface tem- 
perature is reviewed in this report. The turbulent and transitional heat-transfer coeffi- 
cients are based on the approximate solutions of the equations for the thermal boundary 
layer. The laminar coefficients are those due to reference 4 with the corrections for 
variable wall temperature due to reference 5. The boundary layer transition criteria 
used in determining the heat-transfer coefficients are based on the boundary layer mo- 
mentum thickness Reynolds number. Reference 6 discusses data from the literature and 
presents an empirical expression relating the start of transition to the critical momen- 
tum thickness Reynolds number. These Reynolds numbers are presented as a function 
of a gas stream pressure gradient parameter and the gas stream turbulence intensity 

SYMBOLS 

constant 

constant, A^Pr"*^ 
variable used in eq. (5) 
specific heat 

diameter of leading edge of vane 

2 

Euler number, - (dPg/dx)/(pgUg /x) 

Nu/ReJ- ® 

heat-transfer coefficient 
thermal conductivity 
exponents 
Nusselt number 
pressure 


level. 


B 

b 

Cp 

D 

Eu 


■‘^fam 

■-g 

k 

m,n 

Nu 

P 


2 



Pr Prandtl number 

R radius of leading edge 

Re Reynolds number 

SL local Stanton number, h_ /p„C u„ 

X gjj / s p s 

T temperature 

temperature arbitrarily taken from vane temperature distribution 
Tu turbulence level 

Tu mean turbulence level 

Ug local free -stream velocity 

local critical velocity 

X local distance from leading edge of vane 

6* enthalpy thickness 

6 momentum thickness 

X pressure gradient parameter 

p dynamic viscosity 

u kinematic viscosity 

p density 


Subscripts: 

crit 


Heq. 3) 
fam 

m 

rat 

ref 

s 

T. 

1 

t 


critical location of flow transition 
effective gas 
local gas 

gas condition for eq. (3) 
laminar flow conditions 
metal 

correction factor on constant surface temperature coefficient 

reference 

static 

turbine inlet 
turbulent 


w 


wall 


X local condition 

0 refers to positions where flow ceases to be laminar or where it becomes 

fully turbulent 


ANALYSIS 

The design of turbine vanes and blades for gas turbine engines requires accurate 
convective heat-transfer coefficients around the entire airfoil surface. Around the ex- 
terior cylindrical leading edge of the vane, the cylinder correlation is commonly used. 
For the remaining surfaces of the vane, laminar, transition, and turbulent coefficients 
must be supplied. 


Leading Edge Region 


Around the leading edge of a turbine vane the heat-transfer correlation for a cylin- 
der in cross flow (ref. 7) gives 


h^= 1. 14 Re®- ^Pr*^- - 
S D 



0 < |e| < 80° 


( 1 ) 


where Re = pu_D//u and p is determined from the perfect gas law. All properties for. 
s 

the cylinder correlation are evaluated using air properties at the reference temperature 
(ref. 8): 


Tref = 0. 5 Tjjj + 0. 28 Tg + 0. 22 Tgg (2) 

Recent tests have shown that leading edge turbulence could increase the heat -transfer 
coefficient by a factor of up to 1.8 (ref. 9). This effect can be estimated (ref. 10) if the 
turbulence characteristics are known. The correction for this effect, called an augmen- 
tation factor, acts as a factor multiplying equation (1). 


Laminar Flow Region 

For the laminar flow region, the local convective heat-transfer coefficients are ob- 
tained from previous results (ref. 4): 
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(3) 



= F/i 



where 


£am 


Nu 


l/3«^0. 5 


•'Re 


Re 


X 


pUgX 

M 


The “ ^iam ^ evaluated as a function of the local Euler number Eu and 

the local ratio of T^/T^. The values of or Nu/Re®‘^, are approximated by 

solutions for wedge -type flow. These wedge solutions may be used to approximate heat- 
transfer coefficients along an arbitrary profile. The value of Eu at any position along 
the profile determines the corresponding wedge for which, at the same distance x from 
the stagnation point, the heat transfer on the wedge and the arbitrary profile are as- 
sumed to be equal. 

The convective heat-transfer coefficients obtained from equation (3) are based on 
constant wall temperature. Since turbine blades have varying wall temperatures, it is 
necessary to correct the convective coefficients to account for such variations. A meth- 
od for making such a correction exists (ref. 5). The method determines the ratio of the 
variable surface temperature coefficient to the constant surface temperature coefficient: 


\ = ^rat^g 


(4) 


where h^^^. 


is a function of b and b is determined from 



and Tj^ is the local surface temperature, T^ is a reference wall temperature, and x 
is the distance from the stagnation point. One set of data of h^^^. against b has been 
curve fit as follows (ref. 5): 
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(5a) 

(5b) 

(5c) 


= -0. Ill b^+ 0.318 b+ 1.0 0<b^0.6 

h^at = -0- 520 b^ + 0. 310 b + 1. 0 -0. 5 < b < 0 

= -4. 977 b^ - 4. 288 b - 0. 186 -0. 76 ^ b < -0. 5 

Turbulent Flow Region 

The technique for predicting turbulent local convective heat-transfer coefficients is 
based on the approximate integral solution (ref. 3). The solution of the equation for the 
thermal boundary layer takes into approximate account the effects of thermal history, 
free stream velocity variation, and free stream density variation. 

In the solution, the integral energy equation was arranged into a form involving the 
local enthalpy thickness 6*: 

d6* r 1 1 du„ 1 d(T „ - T„)“l h^ 

_Z+ 1 Lw s^ 5* ^ ^ - (g) 

dx [pg dx Ug dx (T^-Tg) dx J pgCpUg 

If a relation could be found between St = h /p_C u_ and 6*, equation (6) could be 

X X D p S f 

solved for 6* as a function of x and therefore h^^ as a function of x. For a zero 
pressure gradient, constant property flow over an isothermal flat plate, the Stanton 
number has the followii^ form: 

= A„Re;:“Pr-" (7) 

The values of A^, m, and n depend on whether the flow is laminar, transitional, or 
turbulent. Substituting equation (7) into equation (6) gives the functional form of the 
Stanton number for the flat plate, 

St = b1/(1-“>(1 - (8) 

-X g* 

where B = A Pr“^. Furthermore, the St in equation (6) was proposed to be some 

VJ A 

general function of Reg* but otherwise independent of body shape, surface temperature 
variation, free-stream velocity, and density variation. If this is true, a universal rela- 
tionship between St^ and RCg* must be given by equation (8). If equation (8) is sub- 
stituted into equation (6) and th? resulting expression integrated, the following result 
will be applicable: 
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Ps^p“s 




( 9 ) 


IB 


l/(l-m)(i _ jj,)l/(m-l) 


/ 


^ p u (T - T )l/(l-m) 
s' s w' 


I m 


m/(m-l) 




where Xq = 

To evaluate B and m in equation (9), the expression for the turbulent local con- 
vective heat-transfer coefficient for a flat plate is used; namely, 

h = 0. 0296 Re°* ^ (10) 

^ * X 

Expressing equation (10) in the form of equation (8) gives B = 0. 0296 Pr"^^^ and 
m = 0.2. Assuming the free-stream dynamic viscosity ju is only a weak function 
of X, 



The factor (T„/T )*^* in the numerator of the right side of equation (11) is the tern- 
perature ratio correction for property variation across the boundary layer (ref. 11) for 
Tw/Ts < 1. The temperature-dependent fluid properties, then, are evaluated at the 
local free-stream static temperature. 


Transition Flow Region 

Over some portions of the surface of a turbine vane, the gas stream flow fluctuates 
intermittently between laminar and turbulent flow. This region of transition is not really 
described accurately by either equation (4) or equation (11). Thus, an empirical ex- 
pression analogous to the form of equation (10) is sought which is applicable in the tran- 
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sition region. Experimental results for heat transfer in the transition region along a 
flat plate have been obtained (ref. 3). For Reynolds numbers Re > 200 000 the heat- 
transfer correlation in the transition region had the form 

Nu^ = 0. 000386 (12) 


Assuming again the form of equation (8), with Pr = 0. 7, the values 
B = 0.000435 and m = -1/9 result. When these values of B and m are sub- 
stituted into equation (9), an expression as a function of x in the transition re- 
gion results; ^ 



To use equations (11) and (13) it is necessary to evaluate 6* at Xq. Since 6* is con- 
tinuous at Xg (between laminar and transitional and between transitional and turbulent) 
it is possible to evaluate 6* by rearranging equation (8): 







B 


1/m 




l-m\psCpUg 



(14) 


At Xq, the values of m, B, and h are from the laminar and transitional expres- 

^XXq 

sions. They are used in equations (13) and (11), respectively. For equation (11), 

B = 0. 000435 Pr'^/^ and m = -1/9. For equation (13), B = 0. 332 Pr"^/^ and m = 0. 5 
when = 0. 332 in equation (3). 

Equations (11) and (13) are finally rewritten in useful form givii^ for the turbulent 
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and for the transitional h 



p^CpU^ (0. 000436) Pr-2/3 (t^ . ^ 

\ w/ 


0. 25 




0.9 


dx + 


1. Ill (2) (0. 332)^ Pr’^/^ ^ 

Sy 


0.000435 




0.9 


-0.111 (16) 


J Xr 


Transition Criteria 

The change from laminar to transitional flow and from transitional to turbulent flow 
occurs when the Reynolds number becomes sufficiently high and thus allows instabilities 
in the boundary layer to grow. To account for local variations in the boundary layer, a 
Reynolds number based on the momentum thickness is used to determine the values of 
Xq in equations (15) and (16) where changes between flow regimes occur. An approxi- 
mate expression for the momentum thickness can be used (ref. 12): 
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(17) 


e = 


0.67 II 


0.5 

s 


0.5 3 



Thus, the Reynolds number based on momentum thickness can be found as a function 
of x: 



(18) 


The values of Xg (where changes occur from laminar to transitional flow) can then be 
found from equation (18). This equation is assumed to apply also in the transitional re- 
gime and up to the location where the flow goes from the transitional to turbulent. The 
question arises as to what to use for Ren • Empirical expressions are given (ref. 6) 

OXg 

for the value of Re^ at the start of transition (that is, the change from laminar to tran- 
sitional flow) for a zero pressure gradient case and a nonzero pressure gradient case. 
They are, respectively. 


Re^ = 190 + exp (6.88 - 1.03 Tu) (19) 

Reg = [0.27 + 0.73 exp (-0.8 Tu)] [550+ 680 ( 1 + Tu - 21X)"^] (20) 

where 

9^ du 

X = ^ 

u dx 

_ rms velocity fluctuation x 100 

and the mean turbulence level Tu characterizes the flow throughout the history of the 
boundary layer. For X = 0, equation (20) agrees with equation (19) within ±5 percent. 
Very little information exists as to what value of Re g to use when the change from 
transition to turbulent flow occurs, but a value of Reg = 360 is used (ref. 11). 
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RESULTS AND DISCUSSION 


A method which calculates laminar, transitional, and turbulent convective heat- 
transfer coefficients over the surfaces of turbine vanes is presented. A sample calcu- 
lation was made for the suction surface of a J-75 size turbine vane whose profile is 
identical to that given in reference 13. The calculation was made for given gas stream 
conditions of turbine inlet temperature equal to 1273 K, turbine inlet pressure of 45. 7 
newtons per square centimeter, fuel/air ratio of 0. 0, vane channel inlet Mach number 
of 0. 31, and a vane channel outlet Mach number of 0. 84. The vane metal temperature 
and free-stream dimensionless velocity distributions used in the calculation are given 
in figure 1. The gas stream conditions given previously were used in the quasi -three- 
dimensional computer program of reference 14 to calculate the velocities. A short com- 
puter program was written to perform the calculation of the heat -transfer coefficients. 
These velocities were input directly into this program. 

The calculated convective heat transfer coefficient distribution for the vane suction 
surface is shown in figure 2. Also shown in the figure are the distributions obtained 
from the turbulent flat plate equation and from a finite -difference computer program 
(ref. 2). The properties of air were used in all cases. 

The flat plate distribution is for a flow assumed to be laminar over the cylindrical 
leading edge and completely turbulent over the remaining vane surface. While closer 
agreement would be obtained by using the transition criteria with the gas stream velocity 
constant, the assumption shows the difference between the approximate integral solution 
and the more conservative (all turbulent) flat plate solution. The approximate integral 
solution and the finite-difference solution both show distinct regions of laminar, transi- 
tional, and turbulent flow. 

The value of the momentum thickness Reynolds number for transition from laminar 
to transitional and from transitional to fully turbulent flow was assumed to be 200 and 
360, respectively, for the integral solution. As a second case of the integral solution, 
the flow was assumed to be fully turbulent near 4. 2 centimeters. This position corre- 
sponds to a film -cooling slot on the turbine vane profile of reference 13. The assump- 
tion of fully turbulent flow was based on the results of reference 15. Experimental tem- 
perature data were presented showing a film cooling layer causing a laminar or transi- 
tional flow to become turbulent. The finite -difference solution requires no a priori as- 
sumption regarding transition. The turbulent viscosity and thermal diffusivity are cal- 
culated from the turbulence model (turbulent kinetic energy and mixing length) and added 
to their respective molecular values. For trq.nsitional boundary layers, the turbulent 
and molecular components are the same order of magnitude. 

For purposes of comparison assume that the finite -difference solution gives the 
baseline heat-transfer coefficient distribution. The turbulent flat plate solution (the con- 
servative solution), except for the vane surface between 4. 0 and 6. 4 centimeters where 
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it is less than 10 percent different, is always 50 to 150 percent higher than the baseline. 
If laminar flow had been assumed at the beginning of the vane, the solution would have 
been lower than the baseline consistently. The approximate integral solution method, 
which accounts for the laminar, transitional, and turbulent flow regimes, produced a 
distribution which is within 15 percent of the baseline between 0. 6 and 2. 2 centimeters, 
where both programs predict that laminar flow exists. Between 2. 2 and 4. 2 centime- 
ters, the integral solution returns transition flow heat -transfer coefficients an average 
50 percent higher than baseline. The integral solution which is fully turbulent at the 
film cooling slot agrees with the baseline within 15 percent from 4. 2 to 6. 4 centimeters. 
The other integral solution does not become fully turbulent until 5. 2 centimeters. From 
that point until 6. 4 centimeters, the two integral solutions are in essential agreement. 
The baseline finite-difference solution goes to fully turbulent flow at 4. 2 centimeters. 
The programmed solution accounted for a film -cooling blowing rate in its calculations. 
The surface region between 0. 5 and 5. 2 centimeters includes the two points of chaise 
(laminar to transitional and transitional to fully turbulent) and the transitional flow re- 
gime itself. The use of the laminar expression (eq. (18)) for the momentum thickness 
Reynolds number in the approximate integral solutions produces an early change to tran- 
sitional flow and a later change to turbulent flow than is produced by the finite -difference 
solution, which uses the calculated boundary layer momentum thickness to arrive at the 
local momentum thickness Reynolds number. The expression for the transitional re- 
gime heat-transfer coefficient (eq. 16), while being of the same level, does not follow 
the finer details of the heat-transfer coefficient produced by the finite -difference solu- 
tion near 1. 0 centimeter. Figure 2 does imply that the heat-transfer rate to the vane is 
about equal for the integral and finite -difference solutions. 

One possible source of error is the assumption that equation (18) is valid in transi- 
tional flow region and up to a value of 360. No assumption at all is made about the accu- 
racy of Reg as calculated by equation (18) once the value exceeds 360. In addition, no 
provision for retransition is made. 

For the variable surface temperature correction of equations (5a) to (5c), if b is 


greater than 0. 6 or less than -0. 76 it is set equal to the respective limit and h^^^^ 


is 


calculated. For the distribution of figure 2, all values of b were within the specified 


range . The correction for approximating wedge flow solutions is made according to 
equation (3). The correction is a function of Eu and T„/T^. If the interpolation for 
^fam the range of tabulated values, is set equal to 0. 332. Equation (3) 

then degenerates to the traditional laminar flat plate solution. In figure 2 this degenera- 


tion occurred at x = 0. 8 centimeter. 


The relative value of the approximate integral solution over the turbulent flat plate 
or the finite -difference solution lies with the user. For the user interested in a conser- 
vative order of magnitude value of heat-transfer coefficient, the turbulent flat plate 
solution and a desk calculator will satisfy his needs. The user of the finite -difference 
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solution program has the greatest flexibility; he will be able to solve the most difficult 
problems, such as the prediction of the local surface coefficient of heat transfer near 
film cooling slots and on the surface of full coverage film cooled vanes. The approxi- 
mate integral solution method gives the user a heat -transfer-coefficient distribution 
which takes into approximate account the laminar, transitional, and turbulent flow re- 
gions. In situations which do not involve blowing into the boundary layer or severe com- 
pressibility effects, the user can obtain heat-transfer-coefficient distributions of 
acceptable accuracy with relatively little effort and with no need to deal explicitly with 
turbulent transport parameters. In addition, when adequate computer facilities are not 
available, the integral technique can be adapted to desk calculators, and it is even fea- 
sible for slide rule calculations. 


SUMMARY OF RESULTS 

An approximate integral solution method for calculating laminar, transitional, and 
turbulent convective heat -transfer coefficients over the surfaces of turbine vanes is de- 
scribed. An example calculation of heat-transfer coefficients using a turbulent flat 
plate solution, the integral solution, and a finite-difference solution was made for a 
turbine vane. The results are as follows: 

(1) The approximate integral solution and the finite -difference solution agree well in 
the laminar and turbulent regions. 

(2) The finite -difference solution predicts a transition region whose length approxi- 
mates that predicted by the integral solution. Its onset begins much sooner, however. 

(3) Use of either the integral solution or the finite -difference solution in a heat flux 
calculation would result in substantially lower levels of heat flux to the vane than would 
the use of the turbulent flat plate solution. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 17, 1977, 

198-10. 
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Vane surface distance (from leading edge), cm 

Figure 1. - Vane metal temperature and dimensionless velocity 
profiles. 
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